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Abstract: The flux of neutrinos originating from cosmic ray interactions with mat- 
ter in the Sun has been calculated based on Monte Carlo models for high energy 
J> \ particle interactions. The resulting flux at the Earth (within the Sun's solid angle) 

is higher than the corresponding one from cosmic ray interactions with the Earth 
. atmosphere. The smallness of the absolute rate, however, precludes it as a practical 

'standard candle' for neutrino telescopes and limits neutrino oscillation searches. On 
the other hand, it facilitates dark matter searches based on neutrinos from neutralino 
■ annihilation in the Sun. 

^5 , 

^ 1 Introduction 



High energy cosmic ray particles, mainly nucleons, interacts with matter to produce secondary 
particles that are mainly mesons and some baryons. This applies in particular to cosmic rays 
entering the atmosphere of the Sun. The produced particles propagate through the Sun until they 
either decay or make secondary interactions producing new particles that contribute further to 
develop a cascade. Decays of particles in such cascades will produce neutrinos and other leptons, 
e.g. muons that in turn decay into neutrinos. This scenario is similar to the cascades induced 
by cosmic rays in the Earth's atmosphere, which we have studied extensively jjj]. However, the 
solar atmosphere is less dense at the typical interaction heights and, therefore, a larger fraction 
of the mesons will decay instead of interacting. This leads to relatively more neutrinos produced 
in the Sun as compared to the Earth. 

It has been proposed ||] that the Sun might be used as a 'standard candle' for neutrino 
telescopes, which is only possible if the flux is significantly higher than the Earth's atmospheric 
flux. A large such neutrino flux from the Sun would, on the other hand, be a severe background 
for searches of neutrinos from neutralino annihilation in the Sun || . The hypothetical neutrali- 
nos appear in theories based on supersymmetry (SUSY) [|J] and are of fundamental interest in 
particle physics as well as in cosmology since they could contribute to the dark matter in the 
Universe. 

In this paper we study the production of muon and electron neutrinos in cosmic ray in- 
teractions in the Sun, as well as their propagation through the Sun and to the Earth where 
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they could be detected in neutrino telescopes, such as Amanda Baikal Dumand 
and Nestor ||. The cascade interactions in the Sun are treated in detail using Monte Carlo 
methods to simulate the high energy particle interactions. In particular, the Lund model Q and 
Monte Carlo programs Jl(J are invoked. The resulting neutrino fluxes at the Earth are compared 
with the fluxes from cosmic ray interactions in the Earth's atmosphere [|l|]. This solar neutrino 
flux is discussed in terms of the above 'standard candle' idea and neutralino search. We also 
investigate the possibility of neutrino oscillations taking place between the source at the Sun 
and the detector at the Earth. 

2 Models and calculational techniques 

2.1 Cosmic ray spectrum 

The flux of primary cosmic ray particles is conventionally parametrised as [l], [Fl], [l2j 

1.7 E- 2 - 7 E < 5 • 10 6 GeV 

(1) 
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The normalisation is here derived [13] from the directly measured primary spectrum using 



balloon-borne emulsion chambers in JACEE |jDf ]. It agrees (within some 10%) with more indi- 
rectly derived spectra based on measured atmospheric muon fluxes |l5f|, and is also compatible 



with the data discussed in ref. [16]. For the energies of our interest (S^lOOGeV), we take the 
cosmic ray flux to be isotropic, since the anisotropy is ^5% j[7]] . 

The cosmic ray composition is dominated by protons with only a smaller component of nuclei 
|l3| . |16[| . Likewise, the outer parts of the solar atmosphere consist mainly of hydrogen with only 
a small fraction of helium. Therefore, the interactions producing the secondary particle fluxes 
are dominantly proton-proton collisions. The small contribution of nuclear collisions can also 
be treated as nucleon-nucleon interactions, since the nuclear binding energies are negligible and 
other nuclear effects have little influence on the high energy secondary particles which we are 
interested in. 



2.2 Solar matter distribution 

The interaction of cosmic ray particles in the Sun can be treated analogously to our calculation 
of cosmic ray interactions in the atmosphere of the Earth [|J. However, the atmosphere of the 
Sun is less dense and particles may therefore propagate deeper in to the Sun such that a more 
elaborate matter density profile is needed. Although there is no well defined solar surface, it is 
convenient to consider different regions in relation to the solar radius Rq = 6.96 x 10 5 km. One 
may then apply an exponential density profile 

p(h)=p e~ h l h \ (2) 

where h > and h < corresponds to locations above and below R & , respectively. The 
parameters po an d ho, given in Table ||, were obtained by fitting the data of ref. [l^] for the 
atmosphere and of ref. [^] for the interior of the Sun. For muons which can reach deep into the 
Sun (~ 20 000 km vertical depth), a third region is introduced to get an adequate description of 
the density profile. 

To calculate the flux of particles from the Sun that reach the Earth one must also consider 
where in the Sun the primary cosmic ray particle interacts. The essential point is whether 
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height h [km] 


po [g/cm A ] h [km] 


h>0 

-2 000 < h< 
h<-2 000 


3.68 • 10- y 115 
3.68 • 10~ 7 622 
45.3 • 10~ 7 2835 



Table 1: Parameters for the solar density profile in Eq. 



the secondary particles and final neutrinos only have to pass through the solar atmosphere or 
through the higher density interior before reaching the Earth. This corresponds to an impact 
parameter b that varies between Rq (peripheral hit) and zero (central hit) and a correspondingly 
varying effective density function along the particle trajectories. An exact treatment of this in 
our Monte Carlo simulation becomes quite complex. We have therefore simplified the situation 
by making the simulations for fixed impact parameter values (b = 0, b = 2R & /3 (geometric 
average) and b = Rq) and then interpolated between them for the integrated results taking 
neutrino attenuation into account. 



2.3 Basic particle interactions 

The interaction of the cosmic ray particles with the solar material is treated as basic proton- 
proton collisions at high cms energy. The production of secondary particles and their decay into 
neutrinos can then be simulated in great detail using the Lund Monte Carlo programs Pythia 
and Jetset |l(J. The neutrino flux arise from the decay of ordinary mesons, mainly tt and K, 
and from decay of muons. 

The production of ordinary hadrons (not containing heavy quarks) is dominantly through 
minimum bias hadron-hadron collisions. The strong interaction mechanism is here of a soft 
non-perturbative nature that cannot be calculated based on proper theory, but must be mod- 
elled. In the successful Lund model || hadron production arise through the fragmentation 
of colour string fields between partons scattered in semi-soft QCD interactions pC(l . The es- 
sentially one-dimensional colour field arising between separated colour charges is described by 
a one-dimensional flux tube whose dynamics is taken as that of a massless relativistic string. 
Quark-antiquark pairs are produced from the energy in the field through a quantum mechan- 
ical tunneling process. The string is thereby broken into smaller pieces with these new colour 
charges as endpoints and, as the process is iterated, hadrons are formed. These obtain limited 
momenta transverse to the string (given by a Gaussian of a few hundred MeV width) but their 
longitudinal momentum may be large as it is given by a probability function in the fraction 
of the available energy- momentum in the string system taken by the hadron. All mesons and 
baryons in the basic multiplets may be produced and the subsequent decays are fully included. 
The iterative and stochastic nature of the process is the basis for the implementation of the 
model in the Jetset program fl~0[| . 

A non- negligible contribution to the inclusive cross section is given by diffr active interactions. 



These are also modeled in Pythia [ 10 1 using cross sections from a well functioning Regge-based 
approach and simulating the diffractively produced final state using an adaptation of the Lund 
string model. These diffractive events are included in our simulations and contribute rather less 
than 10% to the final results. 

In our study of the neutrino flux from cosmic ray interaction in the Earth's atmosphere 
an important point was the production of charmed particles. Their prompt decays gives a non- 
negligible contribution to the neutrino flux at very high energies (E u ^,10 e GeV). This is due 



3 



to the increase in the ir and K decay lengths with increasing energy, such that the probability 
for interacting before they decay increase. The lower density in the solar atmosphere imply 
that their interaction length is larger, such that they still dominantly decay rather than being 
lost through interactions. Therefore, so-called prompt neutrinos from charm decay are not as 
important as for the lepton fluxes in the Earth's atmosphere. 



2.4 Cascade evolution 

To describe the evolution of a cascade in the solar medium we use the same formalism as for 
interactions in the Earth's atmosphere and therefore refer to our earlier study in ref. Q for more 
details. 

The flux of nucleons in the solar atmosphere develops according to the cascade equation 

^ = -^ L + S(NA^NY), (3) 

where X is the depth in the atmosphere and 

^n(E) = P fl (4) 

La °na{E) n A (h) 

is the nucleon interaction length in terms of the number density, n,A,{ti), of nuclei A at height 
h and the nucleon- nucleus inelastic cross section, cjna(E). The first term in Eq. (||) is the loss 
of nucleons due to interactions and S(N A — > NY) is the regeneration due to interactions of 
primary nucleons with higher energies. 

This transport equation is solved through a cascade simulation algorithm as follows. A 
cosmic ray proton is generated with energy drawn from a flat distribution in log E, and a weight 
assigned to reproduce the shape of the primary spectrum Eq. (|]). This primary proton is then 
propagated down through the solar medium according to Eq. (0) without the regeneration term 
S(NA — > NY) resulting in the solution 

<f>{h) = ^o e~ x / XN , (5) 

From this one can obtain the height of the primary interaction by solving the equation 

-\ n ^l = -\ n R= f k dhan(h) (6) 



for h, where R G]0, 1[ is a uniformly distributed random number. 

A proton-proton interaction is then generated in full detail with Pythia [1C] resulting in 
a complete final state of particles. Secondary particles are followed through the atmosphere 
where they decay or interact producing cascades. Secondary nucleons give a flux that is rather 
small compared to the primary flux and could therefore be neglected as a first approximation. 
We do, however, include the main part of this effect by taking into account secondary nucleons 
that have an energy of at least 30% of the primary one. Nucleons with a lower energy give a 
negligible contribution compared to the primary flux due to its steep energy spectrum Eq. (||). 
These leading nucleons emerging in the interactions are therefore allowed to generate a secondary 
interaction at a height calculated according to 

-hx^r = -hxR!= I dhan(h), (7) 
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Figure 1: The E 3 -weighted flux of muon-neutrinos (fn+Vn) and electron-neutrinos (v e -\-v e ) from 
decays of the specified particles (ir,K,[i) produced in cascade interactions in the Sun at different 
impact parameters b. The error bars in (b,e) indicate the statistical precision (bin size 0.1 in 
logiQ E) of the Monte Carlo simulation. Note, the neutrino flux attenuation from propagation 
through the Sun (section 2.6) is not included. 



where H is the production height of the nucleon. The procedure is iterated until the energy 
of the leading nucleon from an interaction falls below 30% of the primary cosmic ray proton 
energy. 

The secondary mesons are propagated down into the Sun until they either decay or interact, 
which is decided by comparing simulated interaction and decay lengths. The interaction length 
is calculated analogous to Eq. (|7|) , while the decay length is given by 

L dec . = -cf3jTlnR. (8) 

Particle decays are fully simulated with daughter particle momenta. In case of interactions, 
the interacting particle is regenerated in the same direction but with degraded energy, chosen 
according to the appropriate leading particle spectrum. Considering only the most energetic 
'leading' particles in secondary interactions is justified because they give the dominant contribu- 
tion to the high energy neutrino flux. Other particles with lower energy are much fewer than the 
particles of the same type and energy produced in primary interactions. Moreover, secondary 
interactions takes place deeper in the atmosphere, where the density is higher, giving a higher 
probability that a produced meson interacts rather than decays. This further decreases the 
importance of particles from secondary interactions. 
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The particle decay-interaction chain is then repeated until all particles have decayed or their 
energy fallen below a minimum of 100 GeV . The energy spectra for neutrinos are finally obtained 
by simply counting the number of generated neutrinos and applying the weight assigned to the 
primary proton. The contributions from decays of the different mesons are shown in Fig. [I] for 
the different impact parameters b. The contribution from charmed and heavier mesons is not 
included, but are unimportant as will be discussed in section ||. 

2.5 Muon propagation in the Sun 

Since muons do not feel the strong force their propagation through matter is quite different 
from the hadrons. Muons interact electromagnetically and rarely experience hard interactions 
with a single large energy loss. Instead they typically loose small amounts of energy in each 
collisions, which can be added such that the energy loss can be treated as a continuous process 
and parameterized in the form 

^ = -ap-P P E, (9) 

with a = 0.0025 GeV/(g/cm 2 ) and (3 = 4.0 • 10~ 6 (g/cm 2 )" 1 . This continuous energy loss 
makes the above treatment of decay versus interaction based on decay and interaction lengths 
unsuitable. Instead, we apply a method where small steps in energy loss (AE/E ~ 10%) are 
taken, under which the energy loss rate is approximately constant (i.e. dE/dh ~ const). The 
decay probability in such a step can then be calculated analytically and Monte Carlo simulated. 
If a muon survives such a step in energy the procedure is repeated until the muon decays or gets 
an energy lower than our 100 GeV cut off. 

The resulting neutrino flux from muon decays are shown in Fig.|l] for the three different 
impact parameters. For muon production with impact parameter b = Rq the effective matter 
thickness is very small resulting in no significant energy loss, such that the muons effectively 
propagate without interactions. 



2.6 Neutrino flux attenuation in the Sun 

The neutrino flux will be attenuated through weak interactions with the nucleons of the solar 
medium. In charged current interactions the neutrino is lost altogether and instead the cor- 
responding charged lepton emerges, which looses significant energy or is absorbed since this 
process mainly occurs in the dense interior of the Sun. In a neutral current interaction the 
neutrino emerges with a reduced energy. We treat this using the analytic technique developed 
for atmospheric cascades of hadrons Q and express the flux by 

- + S{uA^ vY). (10) 



dX A„ 

The first term, the absorption term, gives the loss due to both charged and neutral current 
interactions. The second gives the regeneration due to neutral current interactions of neutrinos 
of higher energies. The latter term is given by 

>^/>'^^^^ 

and can be rewritten as 

Q( A v\ - Mg) [°° ME') K{E) dn uA ^ Y (E',E) _ ^(E) 
S(yA ^ uY )-— ) ] E dE^-^^-^ — -J^E) Z (12) 
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Figure 2: Energy- dependence of the regeneration Z -moment, Eq. (li) for neutrino-nucleon in- 
teractions. 



where the integral defines the regeneration Z- moment [|J. To use this formula it is normally 
assumed that the dependence on the depth X cancels in the ratio c/> v (E')/(j) u (E), which it does 
if the interaction length is only weakly energy dependent. Although this is not the case here, 
this formula is still applicable since the dn/dE spectra are peaked at low energy loss (for high 
energies). Furthermore, one can use the same Z-moment for all neutrinos {v e ,v e ,v^,v^) since 
their spectral form before attenuation in the Sun are approximately the same and their cross 
sections are approximately the same. The Z-moment is calculated with Pythia JHJ using 
a standard technique and it is shown in Fig. ||. Now, Eq. ( [lC| ) can be solved to give the 
attenuated neutrino flux 

Mr) = <t>»fl exp {-(1 - Z)X(r)/\ v {E v ) }, (13) 

where X(r) is the amount of material being traversed. 

The above treatment is of course a simplification. For example, the flux, 4>(E) is a function of 
the impact parameter b. Here, a fit to the flux integrated over the Sun is used. Compared to the 
two most rough simplifications that can be done, neglecting secondary neutrinos and neglecting 
energy loss in neutral current interactions, our simplifications are reasonable. By comparing the 
Z- moments, which are zero if secondary neutrinos are neglected and a NC /o~ NC+cc ~ 2/7 if 
energy loss is neglected in neutral current interactions, we can see that the above treatment is 
a significant improvement. 

In principle, each neutrino should be followed through the Sun with the attenuation folded 
in to get the final attenuated spectrum. However, since the neutrino production only occur in a 
tiny outer fraction of the Sun one may use a simpler factorised approximation. The previously 
obtained fluxes can thus be multiplied by an overall attenuation function for the Sun. This 
attenuation function can then be calculated in terms of the impact parameter 

A(E V , b) = exp {-a(E v ) (1 - Z) X(b)/m N } . (14) 

where X(b) is the effective thickness of the Sun and is the nucleon mass. A numerical 
evaluation of this function is shown in Fig. ^. 

3 Resulting neutrino fluxes 

The unattenuated neutrino flux from the decay of different particles produced in cosmic ray 
collisions in the Sun are shown in Fig. ||. As can be seen from the results for different impact 
parameter values, secondary interactions of mesons and energy loss for the muons are important 
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Figure 3: The attenuation factor, Eq. ^Qj), for neutrinos passing through the Sun as a function 
of neutrino energy and impact parameter b. 



and give substantially less high energy neutrinos at more central interactions. By interpolating 
over different impact parameters and folding with the attenuation factor, Eq. ([HD, the integral 
over the solar disc is carried out resulting in the total neutrino fluxes at the Earth shown (by 
solid curves) in Fig.|3]. The fluxes are compared with the horizontal [21], the vertical [|lj and the 



prompt jlj fluxes from cosmic ray interactions in the Earth's atmosphere. 

As demonstrated in Fig.^J, the fluxes from the Sun are significantly higher than those from 
the Earth's atmosphere. The lower density of the solar atmosphere gives a lower probability 
for secondary interactions and hence favours decays into high energy neutrinos. For the muon 
neutrinos, the solar flux is thus about one (two) orders of magnitude larger than the horizontal 
(vertical) atmospheric flux in the energy range 10 —10 GeV. In the case of electron neutrinos, 
the solar flux very much higher at all energies. This is due to the contribution from muon decays 
(see Fig.|l]), which is only important for the very lowest energies in the atmospheric fluxes. For 
both muon and electron neutrinos, the slope at high energies (<;10 6 GeV) is steeper for the solar 
fluxes compared to atmospheric ones. This is because of neutrino attenuation, Eq. ( |l4|) and 
Fig.[3|, is here getting noticeable in the Sun. 

As mentioned earlier the prompt neutrino contribution, i.e. from hadrons with charm and 
heavier quarks, has not been explicitly studied. The prompt atmospheric fluxes pj are plotted 
in Fig.[| (curve P) for comparison. Due to the short charmed particle life-time, these fluxes only 
depend on the production rate up to ~ 10 7 GeV. The production in the Sun would therefore be 
the same, but the attenuation in the Sun results in a prompt neutrino flux that is lower than 
the atmospheric one. 

Our solar muon neutrino flux is, in Fig.|||, compared with some earlier calculations. The 



results of Seckel, Stanev and Gaisser [22] are in agreement both regarding normalisation and 
shape, although they do not extend to as large energies as our calculation. The result of 
Moskalenko and Karakula [§] shows a significant difference in shape, related to their neglect 
of secondary interactions. By comparing our Fig. [ljab with Fig.[l]c it is obvious that secondary 
interactions and energy loss for muons cannot be neglected, even at relatively low energies. In 
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Figure 4: Cosmic ray induced E s -weighted neutrino fluxes at the Earth integrated over the solid 
angle of the Sun. The fluxes from the Sun obtained in this study (solid lines) are compared with 
the earlier calculation SSG JI?Ij / and the one MK derived from Wj, as well as those from the 
Earth 's atmosphere as calculated for the vertical flux ( curve V) yy, the horizontal flux ( curve 
H) and the prompt charm-induced flux (curve P) 



fact by using the fluxes in Fig. [T]c (where secondary interactions are unimportant) for the whole 
Sun, the result of Moskalenko and Karakula is roughly reproduced. 

One can represent the neutrino fluxes by the simple parameterisation 

f NoE-^/il + AE), E<E , 
HE) = (15) 
[ N( s E-i'- 1 /(l + AE), E > E . 

Although the form is the same as in Eq. (19) of [|l| for the atmospheric fluxes, the physical 
interpretation is not as simple here due to the integration over impact parameter and the inclu- 
sion of the attenuation factor. Still, this form gives a good fit to the total attenuated fluxes in 
Fig. |^ resulting in the parameter values in Table|2|. (Nq is not fitted but given by the continuity 
condition at Eq.) 





N 


7 




A 


E 


i 


% 


»n + ^ 


1.3 • 10 - " 


1.98 


8.5 


• 10" b 


3.0 • 10 d 


2.38 


5.1 • 10" 3 




7.4 • 1(T 6 


2.03 


8.5 


• 10" 6 


1.2 • 10 6 


2.33 


5.0 • 10" 4 



Table 2: Values of parameters in Eq. obtained from fits to the attenuated neutrino fluxes 
given by the solid lines in Fig.^. 

Neutrino telescopes measure neutrino fluxes indirectly through the Cerenkov light emitted 
from the muons and electrons produced in charge current neutrino interactions. The rate is, 
therefore, not directly proportional to the flux, but rather to the flux folded with the probability 
that the neutrino undergoes a charge current interaction and that the produced charged lepton 
reaches the detector, i.e. proportional to its range. Thus, the rate in a water/ice Cerenkov 




og 10 (E th /GeV) 

Figure 5: The integrated rate of muon neutrino events in a neutrino telescope, as given by 
Eq. The event rate for neutrinos from the Sun (S) as compared to the ones based on the 

horizontal (H), vertical (V) and prompt (P) Earth atmospheric fluxes. 



telescope is given by the quantity 

R= [°° dE(j)(E)a u ^{E)^L{E). (16) 



JE th "*p 

Here, p is the density of the medium (p ~ 1 g/cm 3 for water/ice), m p is the proton mass such 
that p/m p is the target number density, and L(E) is the range of the lepton. The lepton energy 
is not the same as the neutrino energy, but for the energies in this study it is a sufficiently good 
approximation. The muon range is given by 

1 (Ep + a/P 



^-r^V^) (17) 

with a, f3 from the muon energy loss formula Eq. (||). The results of a numerical evaluation of 
the rate, Eq. (|T^), are shown in Fig. ||. 

The rate that this flux would induce in neutrino telescopes under construction is very low. 
For example, the version of Amanda currently being deployed will have an effective area of 
about 3 • 10 8 cm 2 and a neutrino energy threshold of about 100 GeV. This would give one event 
per two years of running. 

4 Neutrino oscillations 

If neutrinos are not massless there may be oscillations between the different weak eigenstates 
v e , Ua, v T . This phenomenon has been used to explain deficits of low energy solar neutrinos and 
atmospheric neutrinos in the GeV energy range. A number of particle physics experiments at 
accelerators and nuclear reactors have searched for oscillations. For a recent review of the field, 
sec [p^ j. 

If neutrino oscillations indeed occur, they will affect also the high energy neutrino fluxes 
from the Sun. The situation is, however, very different from the case of low energy solar 
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neutrinos, where the Mikheyev-Smirnov-Wolfenstein (MSW) effect [24] is the dominating source 
of oscillations pq| . The origin of the MSW-effect is the different 'index of refraction' for muon 
and electron neutrinos in matter and it only occurs when certain conditions for neutrino energy 
and nuclear number density are satisfied. These conditions are not fulfilled for the high energy 
neutrinos considered here and, since the fluxes are strongly attenuated when passing through 
the Sun (Fig.||), the effect can be neglected. The oscillations could instead be of importance 
due to the large distance from the Sun as source to the detector at the Earth. 

In our analysis of this, we assume that there are two massive neutrino states that mixes. 
The probabilities that a neutrino emitted as flavour i is of flavour t or £' at detection are given 
by 

P vm {E v ,x) = 1 - sin 2 (2fl) sin 2 (j^j (18) 

P VlVt ,{Eu,x) = sin 2 (2fl) sin 2 , (19) 

where 6 is the neutrino mixing angle, 5m? = \m\ — m 2 \ is the difference of the squared neutrino 
masses and x is the traversed distance, i.e. x = D for the Sun-Earth distance. 

Fig.J^ illustrates the effect of the neutrino oscillation as a function of these two basic oscil- 
lation parameters by showing iso- lines of the indicated per cent changes of the neutrino fluxes, 
i.e. in the ratio 

{^(E u )P u ^(E u ,D) + ^ Ue (E u )P Ue ^(E u ,D)}/^(E u ) (20) 

for the muon neutrino and similarly for the electron neutrino. The results are averages over 
energy bins in log(E u ) centered at E v = 10™ GeV with n = 3, 4. ..8 as indicated. Given the 
starting point with a larger muon neutrino flux than the electron neutrino flux, see Fig.|4], we 
show the 10% and 20% decrease of the muon neutrino flux and the 10%, 20% and 50% increase 
of the electron neutrino flux in Fig. ^ (a) and (b), respectively. 

Two main features can here be observed. First, the lower limit of 5m? increases with increas- 
ing energy. Second, all curves coincide at the highest 5m 2 , where the oscillation length becomes 
much shorter than the Sun-Earth distance such that sin 2 (5m 2 D/4E V ) will be averaged out to 
1/2. 

Fig.J^c illustrates experimental limits and acceptable regions. Region A and B are, together 
with a very small region around sin 2 20 ~ 1 and 5m? ~ 5 • W~ s eV 2 , the regions which could 
explain the solar neutrino problem fl25|| . The region marked C is allowed based on observation 
of atmospheric neutrino fluxes |26|]. D is the region not excluded by the LSND experiment ]|27|]. 



The limit E is set by the Goesgen reactor experiment 28] and F by the accelerator experiment 
L AMPF E734 ]^9j . Since these allowed regions are not overlapping, one can interpret (neglecting 
experimental uncertainties) that as either some result is incorrect or there are more neutrino 
species involved in the oscillation. 

Given the low absolute flux of high energy solar neutrinos, the event rate in currently planned 
neutrino telescopes will be too low for revealing searches for, or studies of, neutrino oscillations. 

5 Conclusions 

We have calculated the high energy muon and electron neutrino fluxes arising from the inter- 
actions of cosmic ray particles with the solar matter. Our resulting muon neutrino flux agrees 



with that obtained by Seckel et al. [22], but our result extends a few orders of magnitude higher 
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Figure 6: Iso-Unes for the indicated percentage (a) increase of the muon neutrino flux and (b) 
decrease of the electron neutrino flux due to v^~y e oscillations as functions of the neutrino 
mixing angle 6 and mass difference 5m 2 = \m± — mil- Results, see Eq. $2(\), averages over 
energy bins in log{E u } centered at E u = W n GeV with the indicated n = 3, 4. ..8. (c) Limits 
and not excluded regions based on different experiments, see text. 



in energy. The muon neutrino flux in ref. |2j is in disagreement with both these results, due to 
an oversimplified model where secondary interactions in the Sun are not taken into account. 

These solar neutrino fluxes are one to two orders of magnitude larger than those from cosmic 
ray interactions in the Earth's atmosphere, when integrated over the solid angle of the Sun as 
seen from the Earth. This opens a possibility to use the solar neutrino flux as a 'standard candle' 
for neutrino telescopes. However, one must here also consider the angular spread introduced 
by the charged current interaction producing the detectable muon relative to the incoming 
neutrino direction. This deflection is typically ~ 10° J 10 GeV/ E^ . In addition, the experimental 
measurement of the muon direction also has a limited resolution. Taking these two effects 
together, the solar disc will typically cover less than ~ 10% of the solid angle that has to be 
integrated over. The solar muon neutrino flux, which stays the same, would then have to be 
compared with a factor ten, or more, increased atmospheric flux such that the two would be of 
the same order of magnitude. Under these conditions, one would only have a factor two increase 
towards the Sun and the use of the Sun as 'standard neutrino candle' does not look so promising. 
One should here also be aware of the very low absolute rate of events in a neutrino telescope 
of the size now under consideration. For example, we estimate the rate of neutrino events with 
E v > 100 GeV to be one per year in a detector covering 6 • 10 4 m 2 . 

We have also investigated the potential for observing neutrino oscillations during the passage 
from the source in the Sun and a detector at the Earth. In principle, one could access an 
interesting region in the parameter plane of sin 2 26 and 5m 2 , but with the very low absolute 
rates it is beyond present neutrino telescopes. 

The positive consequence of these small solar neutrino fluxes is that they will cause less of a 
background problem in attempts to detect neutrinos from other sources. Of particular interest 
here is the search for neutrinos from neutralino annihilation in the Sun, where the predicted rate 
can be up to an order of magnitude larger depending on the supersymmetric parameters [pC[|. 
A clear observation of this phenomenon would both demonstrate supersymmetry, i.e. physics 
beyond the standard model in particle physics, and the presence of non-baryonic dark matter 
in the Universe. 
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